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The effects of electron-electron interactions on tunneling into the bulk of a two dimensional 
electron system are studied near the integer quantum Hall transitions. Taking into account the 
dynamical screening of the interactions in the critical conducting state, we show that the behavior of 
the tunneling density of states (TDOS) is significantly altered at low energies from its noninteracting 
counterpart. For the long-range Coulomb interaction, we demonstrate that the TDOS vanishes 
linearly at the Fermi level according to a quantum Coulomb gap form, viw) — CQ\io\/e'^, with Cq 
a nonuniversal coefficient of a quantum mechanical origin. In the case of short-range or screened 
Coulomb interactions, the TDOS is found to follow a power-law \tjj\" , with a proportional to the 
bare interaction strength. Since short-range interactions are known to be irrelevant perturbations at 
the non-interacting critical point, we predict that upon scaling, the power-law is smeared, leading to 
a finite zero bias TDOS, i'{uj)/i'{0) = 1 -|- (|tj|/ciJo)^ where 7 is a universal exponent determined by 
the scaling dimension of short-ranged interactions. We also consider the case of quasi-lD samples 
with edges, i.e. the long Hall bar geometry, and find that the TDOS becomes dependent of the Hall 
conductance due to an altered boundary condition for diffusion. For short-range interactions, the 
TDOS of a quasi- ID strip with edges is linear near the Fermi level with a slope inversely proportional 
to pxx in the perturbative limit. These results are in qualitative agreement with the findings of bulk 
tunneling experiments. We discuss recent developments in understanding the role played by electron- 
electron interactions at the integer quantum Hall transitions and the implications of these results 
on dynamical scaling of the transition width. We argue that for long-range Coulomb interactions, 
the existence of the quantum Coulomb gap in the quantum critical regime of the transition gives 
rise to the observed dynamical exponent z = 1. 



I. INTRODUCTION 

A. Integer quantum Hall transitions and inadequacies of the noninteracting electron theory 

The physics of disorder and interaction in strong magnetic fields is central to our understanding of the low- 
temperature, quantum mechanical behaviors of novel electronic aiaterials. One of the most important physical 
phenomena under such settings is the quantum Hall effect (QHE).EiH The QHE refers to the low-temperature mag- 
netotrajjisport properties of high mobility two dimensional electron systems (2DES) in a strong transverse magnetic 
ficld.l!rcl The main part of the phenomenology can be summarized by: (1) The existence of stable phases of matter, 
i.e. the quantum Hall states, with vanishing dissipation and integer or fractional quantized Hall conductances; (2) 
The existence of continuous, zero temperature phase transitions between the quantum Hall states, which is often re- 
ferred to as the the quantum Hall (plateau) transitions. The basic physics for the spin-polarized incompressible 
quantum Hall states and their low energy excitations are well understood.au In contrast, (2) is yet an unresolved 
problem, which is the subject of this work. 

In a nutshell, (2) is a metal-insulator transition problem of the Anderson-Mott type in a 2D disordered system with 
strong time-reversal symmetry breaking. These transitions are generally belkived to be prime examples of continuous 
quantum phase transitions, i.e. examples of quantum critical phenomena.ErQ Although there-are reasons to suspect 
that the critical phenomena is universal for both the integer and the fractional transitions J3il2l we focus here on the 
integer quantum Hall transitions (IQHT) in samples with sufficiently strong disorder that fractional quantum Hall 
states do not intervene. In this case, the^transitions are directly between adjacent integer quantized Hall plateaus. 
The experimental data, reviewed in RefJj, can be summarized as follows: (a) On either side of the transition the 
Hall conductivity axy is quantized and the dissipative conductivity has the limit a^x at zero temperature; (b) 
At the transition, axy is unquantized and axx remains finite at zero temperature, so that the disordered quantum 
critical state is conducting. Thus the quantum phase transition is an unusual insulator to insulator transition with 
no intervening metallic phase, only the critical point itself has a finite conductance. 
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In an experimental situation, the divergent length, i.e. the critical singularity is cut off by the presence of a finite 
length scale, giving rise to a finite transition width within which Uxy deviates from the quantized values and a^x is 
nonzero. The transition width, denoted as 5* follows the scaling form, 
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(1.1) 



where L, T and lo are the finite system size, temperature and measurement frequency in a specific experimental 
situation, and Aq, Lq, Tq and loq are microscopic scales. The various exponents in Eq. (1.1) have the usual meaning: 
i/ioc is the static exponent of the single divergent length scale, the localization length ^ ^ o^^°^ where S is the distance 
to the quantum critical point; z^^ is the dynamical exponent defining the length scale introduced by a finite frequency, 
Li^ ~ ^-i/zij. g^jj(j jg |;j^g thermal exponent governing a temperature-dependent length scale Lt ■ 
general, zt and z^^j can b e in dependent exponents,!!^ but zt — fcp-|a|-ffieneric quantum phase transitionJj'B The three 
scaling regimes in Eq. (r]j]_) have all been probed experimentally,t3il3 as well as the regime in which electric field 
strength sets the cut-ofF.tJ The critical exponents extracted from the experiments can be summarized as z/joc — 2.3±0.1, 
^/zuji^ioc = 0.41 ±0.04, and I/zt^^Ioc ~ 0.42±0.04. Thus, we have w/T-scaling with zt = = 1, which is in conformity 
with the dynamical scaling description of a generic quantum phase transition. 

The phase structure of the IQHT appjpjs to be consistent with that of the noninteracting theory of disordered 
2D electrons in a strong magnetic fieldJl3't3 In a single-particle description, all states are localized due to disorder, 
except for those at a single critical energy Ec near the center of each disorder-broadened Landau level. The IQHT 
takes place when the Fermi level i?;- of the 2D electron system and one the discrete set of critical energy Ec cross, 
i.e. when d = \Ep — Ec\ approaches zero. Moreover, numerical calculations based on the noninteracting theory give 
a localization length exponent iy\oc — 2.3, which is remarkably close to the experimentally value. 

However, our understandingpof the IQHT is far from complete. It has become increasingly clear that the noninter- 
acting theory, reviewed in Ref.E3, contradicts in several ways the experimental findings. Three of these are in order, 
(i) Recent experimental work has shown that the tunneling density of states (TDOS) vanishes linearly at the Fermi 
level,Ej in sharp contrast to the finite density of states in the noninteracting theory, (ii) It was pointed out recently 
that due to the peculiar phase structure involving a set of extended states that has measure zero, the conductivity 
axx in the noninteracting theory is rigorously zero in the limit of large sample size at all values of the magnetic field, 
including the critical values, for any nonzero temperaturcEj This is in direct contradiction to the experimental obser- 
vations, (iii) The noninteracting theory does not offer a correct description of the dynamical scaling behavior observed 
experimentally. The dynamical exponent governing how the energy (temperature) scale relates to the length scale for 
noninteracting electrons is z — d — 2^ which disagrees with the experimentally obtained values quoted above. In fact, 
the experimental findings of w/T-scaling with zt = z^j = 1 is in conformity with the dynamical scaling description 
of a gener ic ii qu antum phase transition in which Coulomb interaction is relevant and scales to a finite value at the 
transition.E3Ej 

The failure of the noninteracting theory highlighted by (i)-(iii) puts serious constraints on the ability of the free 
electron model to explain the IQHT in real materials, and necessitates the investigation of the effects of electronic 
interactions and their interplay with disorder and localization. 



B. Recent theoretical developments on the effects of interactions and the focus of this work 

A significant part of the recent theoretical studies on the role of CoulombJintepajCliDns near the IQHT has cen- 
tered around the three inter-connected issues (i)-(iii) raised in the above .oEjSEjEj o There has also been formal 
approaches by Pruisken and collaborators that aim at exteiiding the topological nonlinear a-model description of the 
noninteracting transition to include Coulomb interactions.c3 

The first quantitative study of the one-particle density of state (DOS) in the presence of Coulomb interactions was 
carried out by Yang and MacDonald. Using a self-consistent Hartree-Fock (HE) approach, in which disorder was 
treated exactly while the Coulomb interaction by HE approximation, they found that the TDOS vanishes linearly 
at the Fermi level at all filling factors in the lowest Landau level, even at the critical energyBJ The linear Coulomb 
gap behavior, especially at the critical energy, is in sharp contrast to thp-t expected of the noninteracting theory, see 
(i) above, and is in qualitative agreement with experimental findings.E2l In spite of the dramatic TDOS change due 
to Coulomb interactions, however, Yang, MacDonald, and Huckestein found that the value of the localization length 
exponent and the fractal dimension of the critical eigenstate wave functions remain unchanged from the noninteracting 
theory, as does the qualitative behavior of the conductivity.E3 It is important to emphasize a unique and important 
feature of the HF theory for the IQHT: the noncritical suppression of the single-pacticle DOS, i.e. it vanishes linearly 
at all filling fractions regardless of the whether the system is at criticality or not.c^ 
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In order to understand the effects of Coulomb interactions from the point of view of critical phenomena, Lee and 
Wang carried out a stability analysis of the noninteracting fixed point (NIFP), which governs the noninteracting 
transition, by, numerical calculations of the perturbative renormalization group (RG) scaling dimensions for the 
interactions .E3 They found that interactions of sufficiently short-range are perturbatively irrelevant at the NIFP and 
scale to zero in the asymptotic limit. The NIFP is therefore stable against such screened interactions, and as a 
result, J^ioc — 2.3 and z = 2. Wang, Fisher, Girvin, and Chalkcr have shown that although short-range interactions 
are irrelevant in the RG sense, they generate a nonzero critical value for the dissipative conductanjue, thus remove 
the pathology (ii) of the noninteracting theory, and control the temperature-scaling behavior of axxiLj They showed 
that in the presence of irrelevant interactions, the scaling theory for transport properties becomes u ncon ventional, 



w/T-scaling breaks down and a third independent critical exponent, the thermal exponent zt in Eq. (1.1), emerges. 
The value of zx is set by the scaling dimension —a < of the interaction strength through the finite temperature 
dephasing time in the critical regime, ^ T^p where p — I + 2a/ z, leading to zx = 2z/{z + 2a). They argued that 
quantum critical scaling behavior of this kind may be a generip-|feature of finite temperature transport near quantum 
critical points, when interactions are (dangerously) irrelevant. EJ 

In contrast to short-range, model interactions, true long-range Cpulomb interactions are, on the other hand, found to 
be relevant perturbations at the NIFP, making the latter unstable.E3 Hence the true critical point must be interacting, 
corresponding to an interacting fixed point (IFP) having a finite interaction strength. This is consistent with the fact 
that the experimentally extracted dynamical exponentSj_2jx = z^j = 1, which are typical of the charge dynamics at 
quantum criticality controlled by Coulomb interactions However, Lee and Wang proposed tha|t-the fixed point in 
the theory where Coulomb interaction is treated via the HF approximation may in fact be stable.E2l They introduced 
the concept of a HF fixed point (HFFP) and argued that it is the simplest possible interacting fixed point of the 
IQHT. Correlation effects are found to be marginal perturbations at the HFFP due to the linear Coulomb gap in the 
HF theory that degrades of the RG dimensions of the residual interactions. They conjectured that a change in the 
dynamical exponent (z) with no change in the static one (t'loc) can be due to the noncritical linear suppression of 
the single-particle (tunneling) DOS induced by Coulomb interactions. The HF theory, in particular the HFFP of the 
plateau transition, presents itself as a concrete example. There are two important issues that must be resolved before 
this conjecture can be further substantiated. i— . 

First, the theory of Coulomb gap was derived largely on the basis of classical physics.Ej It applies directly to 
electronic systems with Fermi energy lying in an excitation gap such as semiconductors and insulators. Therefore 
it may not be completely surprising that a 2D Coulomb gap DOS exists away from the transition regime where 
the electronic states are stronghr localized and where the transport is dominated by variable range hopping in the 
presence of a 2D Coulomb gapH. What is remarkably surprising is that the linearly vanishing Coulomb gap is found 
to pertain to the critical regime of the IQHT where the localization length is enormously large and the conductivity 
finite. This behavior is unprecedented and it is nature to ask whether it is an artifact of the HF approximation that 
does not include the screening of the exchange interactions. Therefore, it is necessary to go beyond the HF theory 
in the critical conducting regime and study the behavior of the TDOS when the screening of Coulomb interactions 
is taken into account. In ordinary disordered metals in zero or weak magnetic fields, the dyrpiiucal screening of the 
Coulomb interactions by the diffusive motion of the electrons is known to be very important .Eil^Lj It leads to critical 
corrections of the TDOS.tH The natural question is whether the interplay between quantum diffusion and Coulomb 
interaction at the IQHT leads to a linear Coulomb gap beyond the HF theory. 

This is the focus of the present work. In a recent Letter,Ej we reported our findings that the quantum diffusive 
motion of the electrons, i.e. the diffusive dynamics, is too slow to effectively screen out the Coulomb singularity in 
the dynamical case. A nonperturbative resummation of the the most singular corrections in the long time limit to 
the TDOS gives rise to a linearly vanishing TDOS for the critical conducting state. This behavior, termed as the 
quantum Coulomb gap, can be thought as the quantum mechanical analog of the classical Coulomb gap. It has a 
quantum origin and the slope of the gap is nonuniversal in contrast to the classical case. In this paper, we provide 
more physical and detailed theoretical derivations of the quantum Coulomb gap. We also study the TDOS behaviors 
for short range interactions, both outside the scaling regime where a nonuniversal power-law TDOS is found, and in 
the scaling regime where the power law is smeared and a finite zero bias TDOS recovered in accordance with the 
observation that short range interactions are irrelevant perturbations in the RG sense. Interestingly, in the scaling 
regime, the change in the TDOS, Si>{s) follows a power-law with a universal exponent determined by the scaling 
dimension of short-range interactions and the frequency exponent z — 2 m this case. We also address in this paper 
the issue of whether and how the bulk TDOS depends on the Hall conductance. To this end, we study the case of 
quasi-lD samples with edges, such as in the long Hall bar geometry, and find that the TDOS becomes dependent of 
the Hall conductance due to an altered boundary condition for diffusion in a finite magnetic field. It vanishes linearly 
at the Fermi level with a slope that is inversely proportional to the magnetic field strength in the perturbative regime, 
in good qualitative agreement with recent bulk tunneling cxpcrimcntfO. These results will be summarized in the 
next subsection. 
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The second issue has to do with the impHcations of the hnear Coulomb gap on dynamical scaling. The linearly 
vanishing DOS in 2D means that the averaged energy level spacing scales with the length of system according to 
Ae ~ 1/L, leading to a dynamical scaling exponent z — 1. However, one of the persistent mysteries remains, namely, 
it is not clear that this is the dynamical exponent measured by the transport experiments. The fact that quantum 
diffusion exists at the critical point of the transition implies a frequency-dependent length scale ~ (^^Ticj)"^/^ 
that is shorter than the dynamical length scales derived from the sinpk-particle sector. Notice that the relevant 
DOS in is the thermodynamic DOS or the compressibility dn/d^E3 Although it is somewhat unnecessary to 
associate a critical exponent with diffusion, a value oi z^^ = 2 is directly implied and it should govern the dynamicS|-af 
diffusive transport in the asymptotic limit. In a recent attempt to substantiate our previous conjecture made in Ref.E3, 
Huckestein and Backhaufj evaluated the density-density response function near the IQHT within a time dependent 
HF approximation (TDHFA), in an effort to determine z^^ from two-particle correlation functions. Their analysis gives 
z^ = 1, but under the compressibility sum rule that relates dn/d^ to the static limit of the irreducible density response 
function, it appears to have resulted from using a linearly vanishing dn/dfi in L^. This result is at least counter- 
intuitive, since dn/dfi is expected to be smooth and finite for a disordered system on general grounds. Moreover, 
a finite compressibility is necessary for observing the quantum Hall transition without the latter being interrupted 
by incipient quantization plateaus. If dn/dfj, were indeed vanishing, the linear screening length would diverge and 
the screening, properties of the critical state would be similar to those of an insulator. Recently, Yang, Wang, and 
MacDonaldEZI pointed out that the controversial result may be a consequence of not accounting for the consistency of 
the exchange-local-fields and the disorder potential in the TDHFA used. Analyzing the charge redistribution following 
the insertion of an external test charge, they studied the screening properties in the long wave length limit of the 
self-consistent HF theory and find that the thermodynamic DOS is finite despite of the linearly vanishing tunneling 
DOS in the critical conducting state. Therefore the question of whether or how the vanishing single-particle DOS 
affects the dynamics of transport near the quantum Hall critical point remains open. We will discuss this issue in 
more detail in the last section of the paper. 

The main part of this article is devoted to understanding how linearly- vanishing TDOS in 2D is likely so long as t 
conductivity is finite. Similar analysis in the case of zero naagnetic field has been carried out recently by Kopietz 
in connection to the 2D B — metal-insulator transition.Cj Our basic finding is that in the presence of disorder. 
Coulomb interaction is insufficiently screened by the quantum diffusive medium at finite frequencies. As a result, the 
single-particle DOS in the extended regime comes to resemble that in the localized regime, i.e. exhibiting the linear 
Coulomb gap, although the slope of gap is different due to a different mechanism. After an understanding of the 
dynamics in the single-particle sector has been developed, we will turn to the important question of how the depletion 
of single-particle DOS relates to the larger issue of dynamical scaling near quantum phase transitions in disordered 
systems. 



C. Interplay between disorder and interaction 



At roughly the same time as the discovery of the integer quantum Hall effect, there were some remaiiable devel- 
opments in our understanding of quantum transport such as localization and metal-insulator transitions.ciJ The weak 
localization theory was developed as a perturbative approach to sUady the effects of disorder and interactions. Early 
works by Altshuler and Aronov, and Altshuler, Aronov, and LeeLj found several remarkable effects arising from the 
interplay of interaction and disorder: 1) the electron-electron scattering rate is enhanced due to the prolonged stay 
of electrons near one another; 2) there is a correction to conductivity comparable to the localization effect caused by 
quantum interference; 3) and most dramatically, the TDOS is significantly altered from its noninteracting counter 
part near the Fermi-energy. 

For ordinary disordered metals, perturbative diagrammatic calculations show that in 3D the weak localizatioa. 
correction to the TDOS, I'iuj), is of the form Sv ^ ^JZj, a result largely confirmed by experiments in the early 1980s,Ej 
where uj is measured from the Fermi energy. In 2D, for long-range Coulomb interaction, 

^ - gj^jj In(^To) In(^Ti), (1.2) 

indicating the possibility of a vanishing near the Fermi-energy as w ^ 0. In Eq. ( |1.2[ ), D is the diffusion 

constant, tq is the elastic scattering time, and in terms of the inverse screening length k — 2ne^dn/ d^^ ti is given by 
1/ti — tq(Dk^)^. Summation of all logarithmic terms are needed to find, the limiting behavior. This was done first 
by Finkel'stein in a field-theoretic treatment of disorder and interaction.E^ Defining the dimensionless conductivity in 
units of /h via the Einstein relation a = j^D, it was shown that 
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v{lo) = ,.oe-i;i7 (1.3) 

which is vahd only for a 3> Zn(l/a;ro) such that the weak locaHzation correction to a can be ignored. The conductivity 
is given by cr = vqD where vq is th e finite density of states in the self-consistent Born approximation (SCBA). It 



should be emphasized that Eq. (L3) does not represent the asymptotic behavior of the TDOS at small bias lo where 
the conductivity is strongly renormalized and becomes itself scale-dependent .e3 

The behavior oiv(u)) in metallic systemSpShould be contrasted to the classical Coulomb gap behavior of the TDOS in 
disordered insulators. Efros and ShklovskiiE3 have shown that when the long-range Coulomb interaction is unscreened, 
which is true in dielectric insulators, the single particle DOS exhibits a universal Coulomb gap behavior, 

v^s{oj) = aM'-^le^-^-^l (1.4) 

where d = 3, 2 is the dimensionality and ad is a dimensionless constant of order unity. Thus in the insulating regime, 
one expects a linearly vanishing Coulomb gap in 2D :^es — Oi2\ijj\- Since the long-range 1/r Coulomb singularity is 
crucial in the derivation of the classical Coulomb gap, it is only expected to be valid in the strongly localized regime 
where the screening of the interaction is weak and dielectric like. In the quantum Hall effect, the latter corresponds 

to the regions far away from the quantum Hall transitions. 

Perhaps less well-known is that the same double-logarithmic correction to the XDOS as given in Eq. ( |l.2| ) was 
later derived in the presence of a strong magnetic field by Girvin, Jonson, and Lee£3 and by Houghton, Senna, and 
Ying.LJ Diffusion in ajSirong magnetic field comes from the "skipping" of the semi-classical cyclotron orbits caused 



by impurity scatteringEJ. In SCBA, the diffusion constant in Eq. (1^) is given by Z? = ^''c''o ^' where the cyclotron 
radius Tc = {2N -\- I^^^Ib with £b the magnetic length and N the Landau-lcvcl index. Note that in this case, D 
is proportional to the field dependent scattering rate 1/to{B) ~ [uJc/to{B — 0)]^/^. In this work, we will derive the 



analog of Eq. (L3) in strong magnetic fields by a nonperturbative resummation of the double- log divergences. Since 
the critical conductance is finite and scale-invariant at the IQHT, it is possible for us to derive the true asymptotic 
behavior of the TDOS in the low bias limit. We show that, in the presence of disorder. Coulomb interaction is 
insufficiently screened by the 2D quantum diffusive medium at finite frequencies. As a result, the TDOS exhibits a 
linearly vanishing quantum Coulomb gap behavior. 



D. Main results — quantum Coulomb gap in TDOS 

The main results can be explained physically in a simple semi-classical picture. The electron-electron in- 
teraction v{r — r') can be viewed as being mediated by a fluctuating potential field $(r) with a distribution 
P($) ^ e^5*('")" ^(r-r)*(r)^ ^j^g crudcst level, neglecting all dynamical effects, the presence of a potential 
field directly changes the energy levels of individual electrons and contributes a phase delay e 'io '^^ *(rci(-r )) ^j^^ 
single-electron propagator, where (t') is the classical trajectory. Such a semi-classical phase-approximation has hci 
used recently in the context of composite fcrmions coupled to a fiuctuating gauge field to study edge tunnelingCJt 
Averaging the $-field over different trajectories as well as the random potentials, we obtain the averaged phase lapse 

during a time interval (0,t). This phase delay can be viewed as a Debye- Waller factoiS for the impurity averaged 
single-particle Green's function at a (tunneling) site r, G(r) = (?/'(r, t)V'* (r, 0)), 

G(r)^Go(r)e-^M, (1.6) 

where Gq{t) ~ l/r is the counter-part of G(t) in the absence of interactions. Notice that G(t) no longer depends on 
the coordinate r after impurity averaging. The TDOS is given by 

= -^Im J rfTe'-"^G(T)|,^„_,^+,^. (1.7) 

Two factors, both resulting from the diffusive nature of the electron motion in the presence of disorder, lead to the 
divergence of the phase delay W{t) at large r and subsequently the vanishing of the single-particle DOS at the 
Fermi-energy: 1) electrons stay longer in the vicinity of one another at each encounter, 2) the Coulomb potential is 
not completely screened at finite time. Since we address the effects of interactions only i]p..to. a. phase-delay, this part 



of the physics presumably can be set aside from the rest by performing a C/(l)-rotation.l 

The specific form of the pseudogap in the TDOS depends on the type of the interaction and on the scaling behavior 
of the interaction strength. We shall consider both long-range and short-range, screened Coulomb interactions. 
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1. Long range Coulomb interaction 



In the case of an long-range Coulomb potential, the phase delay diverges at long times as W{t) ^ ln(T/ro) ln(T/Ti). 
We will show in detail that this double-log divergence renders the r-integral over e"'^*^'^^ convergent, whereby enabling 
an expansion in lo for v{uj). In the asymptotic low frequency limit, this leads to a linearly vanishing TDOS at low 
temperatures 

v{uj) = CQh\uj\/e\ (1.8) 

We shall refer to Eq. ( |l.8[) as the 2D quantum Coulomb gap behavior. In contrast to the 2D classical Coulomb gap 
behavior given in Eq. ( p_.4[ ), the coefficient Cq in the quantum Coulomb gap is not a universal number, but rather a 
quantity of quantum mechanical origin. It depends on the microscopic details of the sample au/rh. as the mobility. For 
large w, viiS) crosses over to the perturbative diagrammatic result in strong magnetic fields.BtJ 

Since Coulomb interaction is a relevant perturbation at the NIFP, the true transition must be governed by an 
interacting fixed point (IFP) where the Coulomb interaction strength is finite. Thus we expect that the quantum 
Coulomb gap to be the true asymptotic behavior of TDOS at the integer quantum Hall transitions. It is interesting 
to note that the Coulomb gap TDOS that we |Obtained for the critical conducting state at the IQHT is qualitatively 
different from those obtained in the clean casecZl and in a weak magnetic field. e3 



2. Short range interactions — prescaling regime 



For simplicity, we consider the case of a (5-function interaction potential v{r — r') — u5{r — r') as a prototype 
short range interacting potential. Outside the scaling regime, the scaling behavior of the interaction strength u can 
be ignored, i.e. u can be treated as a constant or equivalently as a marginal perturbation. In this case, dynamical 
screening of the interaction leads to a weaker, single-logarithmic divergence in the phase delay W{t) ~ ln(r/To). The 
integral of the Debye- Waller factor is no longer convergent, such that-a power series expansion in to becomes singular. 
This is similar the situation encountered in the X-ray edge problem.ea We find a pseudogap in the TDOS that takes 
the form of a power-law, 

i.H~t/oKor, (1.9) 

where the exponent a is nonuniversal and depends on the interaction strength. It is well known that transport at 
the quantum Hall transition in the noninteracting theory exhibits anomalous diffusion,^ i.e. the diffusion constant 
D = D{q'^/uj) ^ Do{q'^u))^^'^ when Doq^ > uj, where ry is a critical exponent related to the multifractal dimension 
D2 ^ 2 — 1]. We will show that taking into account the anom alou s diffusion, which has no effect in the Coulomb case, 
only leads to a weak 77-dependence in the exponent a in Eq. ( |l.9| ). 



3. Short range interactions — scaling regime 



Because short range interactions are irrelevant perturbations at the NIFP, the strength of the effective interaction 



u must scale to zero in the scaling regime according to Uch 



UUJ 



x + /z 



, where —x+ ~ —0.64 is the dimension of the 



interaction and z = 2 is the dynamic exponent at the stable NIFP. This makes the phase delay W{t) converge in the 
large-r limit. As a result, the power-law decay in Eq. (1.9) is smeared, resulting in a finite zero-bias TDOS, 



1 



(1.10) 



where luq is a frequency scale and 7 = x^/z ~ 0.32 is a universal exponent. This result leads to several interesting 
predictions: (a) For short-ranged interactions, the TDOS is finite and nonuniversal at zero bias, (b) It can be shown 
that i>{0) <C vq if the bare interaction strength is strong, so short-range interactions irrelevant in the RG sense can 
still lead to strong density of states suppression at low bias, (c) Eq. (I.IC) shows that the TDOS increases with 
LO according to a universal power law with an initial cusp singularity for our value of 7. These predictions can, in 
principle, be tested experimentally since the Coulomb interaction can be made short range by placing a metallic 
screening gate (ground plane) nearby. 
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4- Quasi-ID samples with edges 



We als o study wh ether and under what condition the bulk TDOS depends on the Hall conductance. That 
in Eqs. (1.8, r9| , ni0| ) does not depend on Oxy is a direct consequence of the fact that transverse force does not 
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affect the charge spreading in the bulk of the sample. Thus any direct cra;j,-deperidence in tunneling must come fr^ 
contributions at the boundary. It has been shown by Khmel'nitskii and Yosefin,Eil and by Xiong, Read and Stoni 
that in the presence of edges, the Hall conductance enters measurable quantities even in the perturbative limit. More 
recently, Shytov, Levitov, and Halperin studied the problem of edge tunneling into the fractional quantum Hall state 
where the HalL-conductivity dependence of the I-V characteristics also arises from the boundary condition at the 
tunnehng edgeO 

We considered a quasi-lD sample with its length L much greater than its width W and with two reflecting edges. 
This condition is realized experimentally in the long Hall bar geometry. Edge effect becomes prominent in such limit 
since the boundary condition effectively changes the diffusion constant from D in 2D to Did = D{1 + ^fj) where 
= is the Hall ratio. For the case of (5-function interaction, we find that the asymptotic TDOS becomes linear 
in an infinite strip with edges: 

u{uj) = s\u;\. (1.11) 

The slope of the linear gap is proportional to the inverse of the dissipative resistivity s ~ p~^. In the high field 
limit where energy levels form Landau bands (uIcTq ^ 1), the bare value of a^x in SCBA at the center of the Landau 
levels is proportional to the Landau index A'" while the bare Hall ratio is of order 1. The above result then implies 
that, in the perturbative regime where the localization correction to the conductivity is much smaller than the SCBA 
conductivity, the slope of linear density of states is proportioiial to A'^ or 1/B. Interestingly such a dependence was 
indeed observed in the high field bulk tunneling experiments.Ell It remains to be seen whether the samples used in 
certain experimental setup can be qualifies as being quasi-lD with edges and the presence of ground planes indeed 
make the Coulomb interaction short-range. 



E. Organization of the paper 



.e.role of electronic interactions at the NIFP of the quantum Hall transition, and cast 
on the RG dimensions of interactions in an the analysis of the level spacing in the 



In Section II, we revisit tl 
the results of Lee and Wang 

Hartrecj-Bock theory. In doing so, we make connections to the more conventional scaling theory of localization a la 
WegnerE^I and motivate the study of the single particle DOS. We then proceed, in Section III, to formulate the effective 
field theory and the semi-classical phase approximation, and to derive the Debye- Waller factor in Eq. (1.6). The bulk 



TDOS in 2D is calculated in Section IV for various forms of interactions in the perturbative and the scaling regimes. 

,(1.9),(1.1C) are derived in this section. The effect of the anomalous diffusion at 



1. 



The results summarized in Eqs 
the IQHT is also s tudied. Section V is devoted to the derivation of the TDOS in quasi-lD samples in the presence 
of edges, Eq. ( 1.11 ). Summary and discussions on the connection between the single-particle DOS and the dynamical 
scaling of the transition width are given in Section VI. We argue that the existence of the quantum Coulomb gap in 
the quantum critical regime of the transition gives rise, through the interplay between quasiparticle decay rate and 
level spacing, to the experimentally observed dynamical exponent z = 1. 



II. ROLE OF INTERACTIONS — HARTREE-FOCK ANALYSIS OF LEVEL SPACING 



The Hamiltonian of interest describes interacting electrons moving in a two-dimensional random potential in the 
presence of a magnetic field: 



H = 



^ 2m V 



-A 

c 



(2.1) 



where A is the external vector potential producing the strong transverse magnetic field, V{r) is the one-body impurity 
potential. The interacting part of the Hamiltonian i?i„t is given by is the two-body interaction potential. 



U 



|r, -r,|f 



(2.2) 
(2.3) 
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Here u and p control the strength and the range of the interaction. The Coulomb potential corresponds to m = e^/e 
and p = 1, and a short distance cutoff can be introduced for the case of p > 4. 

The noninteracting fixed point of the IQHT is obtained by setting u — 0. The effects of interactions can be studied 
in the framework of critical phenomena by analyzing the stability of the NIFP. Imagine starting with a system at the 
NIFP, and adiabatically switching on the interaction u. One can ask whether u is a relevant or irrelevant perturbation 
in the renocmalization group (RG) sense by calculating the RG scaling dimension of u. This has been done by Lee 
and Wang.E3 They analyzed the most singular contributions to the disordered averaged free energy functional. Here 
we present an alternative treatmento from the point of view of the single-particle DOS and cast the result of Lee and 



Clfi- I 

chE 



evel spacing. 

at criticality. Consider a finite system of 



Wang in terms of the interaction induced corrections to the single-particl| 

For this purpose, it is convenient to use the exact eigenstates approae 
linear dimension L, and two adjacent one-electron eigenstates located on the two sides of the critical energy with 
energy Ei < Ec and E2 > Ec and a separation oj — Ei — E2. In the noninteracting theory, the finite size scaling 
behavior of the separation should follow that of the mean level spacing and scale with L according to 

where is the noncritical DOS in the noninteracting theory. Switching on the interaction u causes a mixing of the 
critical eigenstates, which results in shifting of the single-particle energy levels. The level spacing becomes, 

A12 = A?2+<5Ai2, (2.5) 

where (5Ai2 is the level shift of E2 due to the mixing with level Ei. The question we would like to ask is, in 
the limit lo ^ and the associated length scale = {vqujY/'^ — > 00, how the interaction correction to the level 
spacing <5Ai2 scales with L as we approach the thermodynamic limit L < — > 00. If it falls off faster than the 
mean level spacing ~ of the noninteracting system, the level statistics will be determined by that of the 

noninteracting eigenstates and unaffected by the interaction in the thermodynamic limit. Thus the interaction would 
be an irrelevant perturbation in the RG sense. On the other hand, if 6A12 falls off slower than the level spacing 

will be dominated by the interaction induced level shifts as an increasing number of noninteracting eigenstates gets 
mixed by the interaction with increasing system size L. The interaction would therefore be a relevant perturbation 
in this case, and the noninteracting fixed point would be unstable. 

Let's calculate 5Ai2 perturbatively, which is sufhcient for the perturbative RG analysis. To first order in perturba- 
tion theory, we have 

JA12 -Pi(2|iJint|2), (2.6) 

where Pi is a projection operator that keeps exclusively the contribution from state |1) to the level shift of |2). The 
factorized interaction iJi„t can be written as 

7Jint = E^"^C"C„. (2.7) 

a 

where is precisely the Hartree-Fock self energy correction to the single particle state |a) — c^|0), 

occ. » 

= E / [\Mr)\'\Mr')\' - MrrMr')mr')Mr)] «(r - r')- (2.8) 

Here the summation is over all occupied states. To discuss the average energy shift, it is necessary to study the 
disorder average of the self-energy at the fixed energy E2, 

^2 = ^E^(^2-£;a)Sa, (2.9) 



where vq is the DOS per unit area for noninteracting electrons. Taking the disorder average of Eq. (2.6) and using 



Eqs. (2.7), (2.?), and (|2.9|), we obtain the averaged level shifts. 



SA12 = (— ^) / dWr'^ [\MrW\Mr'W - Mr)*Mr')r0{r')Mr)] 5{Ei - E2 ~ oj)vlv ~ v'). 

= J d^^d^^' [OA{r - r') -Oeir- r')] i^(r - r'). (2.10) 
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In order calculate this quantity, we need to kn ow the scaling behavior of the impurity average products of four wave 



functions denoted by Oa and Ob in Eq. ( 2.1C ) in the limit w ^ 0. The latter are functions of r — r' since translation 
symmetry is restored after impurity averaging. ___ 

It is instructive to follow Wegner's symmetry decompositionEj^E^EZl and extract the orthogonal (cigcn) scaling 
variables under the RG. Consider the most general four-field operator in the unitary ensemble, 

o,= y: ^11 ^ E (2-11) 

where O^^ — O^^. The coefficients v obeys the traceless condition, corresponding to the subtraction of vacuum 
expectations, 

E«:? = E-^] = o. (2.12) 

There are two irreducible representations for these operators, a symmetric one and an anti symmetric one, under 
permutations of indices of the rank two tenser. We can therefore decompose O4 into independent scaling operators, 

O^J, = \o+ + \o-, 0±=Oi\±0'J'^. (2.13) 

The operators 0+ and 0_, having independent scaling dimensions x±, describe the eigen scaling directions of the 
four-field operators under the RG. In the unitary universality class of the metal-insulator transition i.e. the cases of 
weak magnetic field and spixL-flip scattering by magnetic impurities, x± = ±^/2e has been derived from perturbation 
theory in 2 -I- e dimensions .Lj r~ir~i 

At the IQHT, the RG dimensions of 0± have been determined numerically by Lee and Wangj£3'CJ from the leading 
scaling operators associated with the fusion products of four fermion operators that are antisymmetric and symmetcie 
under permutations respectively. The scaling dimension of 0_ is obtained from the product of two "spin" operatorsE3 

x_ = X2s = -0.60 ± .02, (2.14) 

whereas tkat of 0+ is extracted from the leading scaling operator fused by the product of two nearby density 
operators^ 

a;+ = a;2p = 0.65 ± .04. (2.15) 



g scaling behavior will be 
is the ensemble averaged 



A general four- field operator involves contributions from both 0_ and 0+, but its leai 
dominated by that of 0_, since 0_ is much morpirelevant than 0_|_. A good exampL 

inverse participation ratio introduced by WegneiH P'^'^\ Expressed in terms of a four-field operator, its scaling 
dimension, which is also known as the multifractal dimension of the eigenstates, £'(2), is governed by a;_, i.e. D{2) = 
d-\-X- = 1.4 ±0.02, indicating strong amplitude fluctuations of the critical eigenstates at plateau transitions. Indeed, 
the exponent rj used by Chalker and DanielE2l to describe the anomalous diffusion at IQHT is given by 77 = —x-. The 
scatter in the value of a;_ is most likely due to uncertainties involved in different numerical approaches. By the same 
token, one can show that the scaling behavior o 



the ultrasonic attenuation, extensively studied in 2 -f e dimensions 

near the conventional metal-insulator transitiond£j, is controlled by 0_ as welL 

Now let's apply these results to the interaction induced level shifts in Eq. ( 2.1(]| ). We will show that the density- 



density correlation is, in contrast to the inverse participation ratio and the anomalous diffusion coefficient, controlled by 
the symmetric operator 0-|_ with the scaling dimension a;+. Individually, operators Oa and Ob contain contributions 
from both 0_ and 0+ . Their leading scaling behavior is therefore dictated by that of the operator 0_ . We have 

0^(r-r')~OB(r-r')^ (^^^y , (2.16) 

for \r — r' I L ^ L^^. However, the combination of Oa ^ Ob has precisely the symmetry of the symmetric operator 
0+. Thus, 

Oa{v - r') - Ob{v - r') ^ {^^^^ ^ ■ (2.17) 
We see that the scaling behavior of the interacting-induced level shift is determined by the fusio n pro duct o f two 



density operators in the symmetric representation. The absence of pure powers of |r — r'| in Eqs. (2.16) and (2.17) 
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comes from the fact that bi-Hnear field operators have dimensions zero,^t£3 consistent with the single-particle DOS 
being non-critical at the noninterac ting critical point. 

Substituting Eq. (2.17) into Eq. ( ^.10| ) and carry out the spatial integrals from a lattice cutoff a to the system size 
L, we obtain, 



(5A 



12 



U 

Up 



(2.18) 



where ci and C2 are nonuniversal constants. Notice that for p > 2 + a;+, the integral depends on the lower cutoff and 
the second term in Eq. (2.18) diverges as a ^ which must be absorbed into the renormalized interaction. We are 



now ready to determine the relevance of the interactions by comparing the scale dependence of (5Ai2 to the mean 
level spacing of the noninteracting system for large L. Define the scaling dimension of the interaction u according to, 



(2.19) 



dlnL 




we obtain using Eq. ( 2.18 ) and Aj'2 ^ L 



max(2 — p, — x+). 



(2.20) 




irrelevant 



FIG. 1. The scaling dimension x of the interaction 
strength u in a l/r*'-potential as a function of p. The 
interaction is relevant for a; > and irrelevant for x < 
in the RG sense. 

The behavior of x is shown in Fig. 1. For p < 2, the interactions acquire a RG scaling dimension x = 2 — p > Q 
and are relevant. We can refer to these type of interactions as long range interactions. cj In this case, the interaction 
induced level shift becomes much larger than the mean level spacing of the critical eigen states in the noninteracting 
theory for large system sizes L. For the Coulomb interaction, p = 1, u has a RG scaling dimension xcoui. = 1 and is 
therefore a relevant perturbation. The resulting flow away from the NIFP will lead to an interacting fixed point at 
which the effective interaction strength is finite. Presumably, the simplest version of the latter is the Hartree-Fock 
fixed point discussed in the previous section. At the level of the Hartree-Fock theory, Eq. ( 2.18| ) shows that the level 
spacing is entirely dominated by the Coulomb interaction induced level shift that scales as in the thermodynamic 
limit, consistent with the linear Coulomb gap DOS found in numerical calculations .cJ 

On the other hand, for all values of p > 2, we have a; < 0. The interactions are irrelevant and can be refereed to as 
short range interactions. The dipole-dipole interaction, in particular, having p = 3, belongs to this class of interactions. 
The NIFP is therefore stable against short range interactions. For screened Coulomb interaction with 2 < p < 2 -I- a;+ , 
the RG dimension of u is x = 2 — p, while for p > 2 + x^, x — —x^. In both cases the interaction scales to zerii-at the 
transition in the asymptotic limit, although, it controls the finite temperature behavior of the conductances. lj From 
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Eq. ( p.lq ), it is clear that the interaction induced mixing between the critical eigenstates of the noninteracting theory 
only leads to level shifts that are much smaller than the mean level spacing in the thermodynamic limit. Thus, the 
zero bias DOS must be finite in the asymptotic limit for short range interactions. However, as we will demonstrate 
later in this paper, even in this case, the interactions lead to remarkable properties of the TDOS in the pre-scaling 
regime, which may have important experimental consequences. 



III. FIELD THEORY FRAMEWORK AND SEMI-CLASSICAL APPROXIMATION 

In order to include the screening of the Coulomb interactions by the diffusive electrons in the calculation of the 
TDOS beyond the Hartree-Fock theory, we will set up, in this section, the effective field theory and the semiclassical 
approximation in order to derive the Debye- Waller factor in Eq. (pT^). 



A. General formalism 



We consider the actio n for interacting electrons in a random potential and a magnetic field described by the 
Hamiltonian in Eq. (2.1), 



drd^vC, 



where 



C = ^* [a^ + Ho + V^(r)] V + ^ j rf^rV'*(r)V'(r)v(r, r')V'* (r')V'(r')- 



(3.1) 



(3.2) 



Here, once again V{y) is the random external potential, u(r — r') is the electron-electron interaction, and TJq = 
2^(9i + ieAjY wi th A(r) being the vector potential of a uniform external magnetic field perpendicular to the 2D 



plane. In Eq. ( |3.2| ), ip and ip* are independent Grassmann fields. The electron single-particle Green's function is 
defined as: 



where Z is the partition function expressed in terms of the imaginary time path integral. 



D[ii;*]D[ilj]e- 



(3.3) 



(3.4) 



The interaction between the electrons can be viewed as being mediated by a fluctuating scalar field $ in the Coulomb 
gauge, for one can always rewrite the interaction term in the following way: 

g-^i / cPr J d2r^*(r)V(r)t>(r,r')V'*(r')V(r') ^ J J^[q,y J d^r^.' W^WV^M- i / d^r J rf=^r' *(r)u- ^ (r-r')*(r') ^ (-35-) 

To perform the averaging over random potential, we use the replica trick, calculate 2", where n is the number of 
replicas, and take the limit n — *■ at the end. The ensemble averaged single-particle Green's function can be obtained 
according to, 



(G(t)) - limy Dm J D[V]P[V] J i?[^*]2?[7/;]V^„,(r,0)C,(r,T) 

C ^.p " r /• 

exp{~ rfrV d^rrair,T)[dr+H„ + V-i^a{r,T)]i:^{r,T) 

+ J dr J dr'<^o.ir,T)v^\r-r')<t>a{r',T) |. (3.6) 

In the above equation, P[V] is the distribution of the random potential which is taken to be Gaussian white noise: 
P[V] ~ e^^^ '•'"^ for the short-range correlated impurities, a is the replica index, and ai represents an arbitrary replica 
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channel. As in the usual treatment of disordered systems, integrating over V{r) in Eq. ( |3.q ) introduces a 4-point inter- 
action term in the action that is local in space but nonlocal in time, fSaSa'/ '^^^ So '^'^ jo '^'^'\'^a{Y-,T)\'^\'ilja'{Y,T')\'^ . 
Th latter is usually decoupled by introducing an auxiliary field Qa,a' {t, t') by way of the Hubbard-Stratonovic trans- 
formation: 

/ d^v /^''dr'|V„(r,r)|^|^„,(r,r')|^ 

D[Q]e-^-l' /o^ fo 'ir'Q^^dr,r')Q^'Ar',r)+^ J d=r J^" dr dr' K (-)Qo„' (-,r')^v (r') ^ ^37) 

This quench- averaging process leads to the following replicated action. 



Sir 



,^,Q,<^)^ f dh f drl f dT'V{C(r,r){[9,+i?o-<i>a(T)]<5(T-T')5o 
J Jo Uo 

- iQa,a'{T,T')}^/ja'{r,T') + —Qaa' (t, T')Qa' air' t)} 



+ / d'r'<i>„(r,rKi(r,r')$„(r',T) 

The impurity-averaged Green's function in Eq. (|3.6|) becomes. 



(3.8) 



(G(r)) = lim y Dm J D[Q] J D[r]D[^]rdr,0)rA^,T)e-^^^''^''^'''\ (3.9) 

The rest of this section proceeds as follows: 1) We integrate out the Q-field as well as the V'-field and derive an 
effective action in terms of the ^-field, Ses{^). The Green's function (G(r, r';r)) can be expressed in terms of the 
averaged non- interacting electron Green's function in the presence of the fluctuating potential which we denote 
by G(<i>), weighted by e"'^"^'*^*^ The effective action ScS can be obtained systematically in a power series of <& and 
I/gxx- 2) By appealing to the semi-classical approximation for the slowly- varying part of the $-field, we argue that 
G($) ~ G(0) exp(— J dr J d'^r^p), where p(r) is the diffusion propagator. 3) Keeping in Scs up to quadratic terms 
in $^ and integrating out the <& field, we arrive at (G) ^ G{0)e~^'^'^\ 

B. The effective action 

Let us define the effective action by formally integrate out the Q and i/'-field, 

Z[$] = e-^=«(*) = J D[Q] J D[r]D[^]e'^^''''''^''^''^\ (3.10) 
Carrying out explicitly the functional integral over the -f/^-field, we obtain 

Z[$] = e-^--(*) = J i?[Q]e^««W-^2"'(*), (3.11) 

where, 

Sgi"^) = -:^Tr02 + Trln[icj ~ Ho + i<i> + iQ], (3.12) 

4'')($) = iTr$(r)w-i(r,r')$(r'). (3.13) 



In Eq. (3.12), uj is the fcrmion Matsubara frequency and "Tr" stands for the trace over the replica, spatial, and 



Matsubara indices. Expanding Trln[iLj — Hq + i$ + iQ] in power series of we get 

Sq = -^TrQ'^ + Tr \n[iuj - Hq + iQ] + Tr (- -J -i$ 



2 \tuj ~ Hq + iQ iuj — Hq + iQ 
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We group the terms in Sq{^) which are zeroth order in ^ as Sa, 



Sa = -T^TrQ' + Trln[iuj - Ho + iQ], 



the terms Unear in <i> as AS'i($) 



AS*! = Tr 



1 



iuj — Hq + iQ 



and the <i>^-term as AS'2($), 



A5, = — Tr 



1 



1 



iu! — Hq + iQ iu! — Hq + iQ 



(3.15) 



(3.16) 



(3.17) 



Clearly, So- is nothing but the transport action for electrons in a random potential in the absence of the Coulomb 
interaction. The standard procedurecj is to expand around the saddle-point of Sa- The self-consistent saddle-point 
equation is given by: 



iQ{r) = -g{r\ 



iuj — Hq + iQ 



The saddle-point solution is given by 



^[Qo]^ 



(3.18) 



(3.19) 



where tq is the elastic scattering time. In general tq depends on the magnetic field B. In the weak magnetic field limit, 
the Landau levels overlap due to disorder broadening and to{B ^ 0) = l/27ri^o3 where vq is the density of states. 
In the strong field limit, uJcTo{0) 3> 1, where Uc is the cyclotron frequency. The diffusion co mes f rom the "skipping" 
of the semi-classical cyclotron orbits caused by impurity scattering. One must solve Eq. (3T8)|i£i the presence of 
quantized Landau levels. Such a solution renders the famous semi-circle density of states by Andcc3, usually referred 
to as the result of the SCBA. Let us denote the saddle-point Green's function as 



[G. 



1 



splaf) 



Hq + iQo 



|r)<5Q/3(5mn, 



(3.20) 



which is also called the SCBA Green's function. Using Ggp one can calculate the bare parameters of the theory, i.e. 
the transport coefficients axx and[«iWi^f^^ level of SCBA. Without going into the details of the SCBA calculations, 
which can be found in referencescSEj^tSH, we point out the following key results: 1) in the strong field limit the 
elastic scattering time tq is of order ■\/to(0)/cl'c, the mean- free path £ = w/tq becomes the cyclotron radius Rc and 
the diffusion constant D — ^R'^/tq depends on the magnetic field; 2) axx in the center of the A^th Landau level is 
approximately given by cr^x^^ — [N -f l/2)e^//i, therefore the perturbative expansion in \/axx is valid as long as 
> 1. 



The single-particle Green's function in Eq. ( 3.18 ), and hence the saddle-point solution for Q has, quite generally, 
a branch cut at w ~ 0. Taking this into account, the saddle point solution assumes the following generic form, 
iQo = 90 + ^0^1 where A is a diagonal matrix, A^™ — 5„m(5Q,/3sign(r7,), in the space spanned by the replica (a, (3) and 
the Matsubara frequency (n, m) indices. The non-linear a model is obtained by including the gapless, long-distance 
fluctuations around the saddle point manifold of the form 



Q(r) -C/-i(r)QoC/(r), 



(3.21) 



where U \s a. unitary matrix U G U{M), and M is the product of the number of replicas and that of the frequencies. 
Defining a new dimensionless matrix field 



Q(r) = f/-i(r)A[/(r), 



(3.22) 



it is straightforward to show that 



A^i = Tn^oTriQ^) + i-Tr[(G,p + Gl^m + 0{QVQ), 



A52 - 7ri.o^Tr($^ ~ Q$Q$) + gTr[(G,p 



G*)<^{Gsp + G*)^+0{Q^Q) 



(3.23) 
(3.24) 
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It is now c onvenient to denote the quantum average over Q under the statistical weig ht e-'^-('3) by (. . We obtain 
from Eqs. ( |3.1l| - ^T7 ) the effective action, 



5eff(<i>) - (A5i). + - \ [{{AS,f). (A5i)^] + 5f ($) + (3.25) 

The second and third <i?^ terms in Eq. ( 3.25| ) define the (density) polarization function 11, 



- i [{I\Sl). ~ (A^i)2] fd'r f dr'n(r,r';^„)*n(r)$-„(r') (3.26) 

The polarization function can be calculated in power series of l/axx- To leading order in l/cr^x, we recover the result 
of the ladder-approximation: 



n(q,ic^„) = ^o ^ 21.. | . (3.27) 



Higher order interference corrections presumably renormalizes the diffusion constant I?, and the thermodynamic DOS 



i/Q dn/dfi in Eq. (3.27). Thus, we have derived the effective action to order <i>^. 



SeS-^Y. W [^"'(^ - + - r'; iuJu)] $„(r'). (3.28) 

n 

The scalar field <& precisely mediates the diffusion-screened electron-electron Coulomb interaction.0 

C. Semi-classical phase-approximation 



Now we turn to the evaluation of theimpiurity averaged single particle Green's function given in Eq. (3.9). By a 
simple reordering of functional integrals,oE3 we have 

(G(r,r,r))= J D[<i>]G{r,r,T,^)e-^''"'-'^\ (3.29) 

where G is the averaged Green's function in a fixed configuration of the scalar potential $, 

G(r,r,$) = hni y D[Q] J Di^D^p^^.^a^e-''^^' ^^^'^^''^ / Z[^]. (3.30) 

We now make an approximation regarding G(<i>) that takes into account exclusively the important interference effects 
between the phases of the electron wave functions. The amplitude fluctuations are small for the slowly varying 



fluctuations of the $-fleld that dominate the contributions to the effective action in Eq. (3.2? ). Since these fluctuations 
are spatially smooth on the scale of the elastic mean-free path £, i.e. V^t/Ef ^ 1, they do not significantly alter 
the classical trajectory of the electrons. This is a unique feature of the slow diffusive dynamics of the electrons in a 
random media. Appealing now to the semi-classical approximation, the single-electron propagator in the presence of 
interactions is modified by a U(l) phase factor: 

G(r,r;T;$) « G(r,r;T;0)e-^^='(*^ (3.31) 
where AS'c; is the change of the action caused by sampling the potential <i> along the classical path: 

AS,m = -I r dT'$[r^,(T'), r'], (3.32) 
Jo 

with r^j(T') being the classical trajectory that starts and ends at r in the presence of random potential but in the 
absence of the $-field. Upon averaging over the random potentials, the classical trajectory can be described by a 
random walk. Let p(r, t) be the probability of a particle being at r at time r, 

'f{Kiir'),r') = J drV(r',T')$(r',r'). (3.33) 
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Since the critical conductivity is finite at the IQHT, the charge spreading is expected to be described by (anomalous) 
diffusion. The probability density p then satisfies the following diffusion equation, 

[-D{Wr + dr'W, r') = [<5(r') - S{t' - r)]<5(r' - r), (3.34) 

where the (5-functions on the right hand side result from the boundary conditions imposed on the original trajectory 
and correspond to injecting an electron at r and time and removing it at time r. The associated current density is 
given by 



-D{V-jHZ X V)p, 



(3.35) 



where — cr^y/axx is the Hall ratio. Notice that does not enter the diffusion equation (3.34) because the 
transverse force does not affect the charge spreading which is described by V • J in the continuity equation. Solving 
Eq. (3.34) in the bulk of systems without edges, we find. 



p{q, iujn) 



(3.36) 



Later we will show that does enter in the presence physical edges. In that case, the diffusion equation must be 
solved with the appropriate sp atial b o unda ry con dition s. 

Inserting the results of Eqs. ( |3.36D , ( |3.32| ) and ( ^.33]) into Eq. ( ^.3l| ), we have 

G(r, r;r;$) w G(r,r;T;0)e Jo J \ . )p\ , j ^ 



(3.37) 



Note 
fields 



±he above is but a special case of the more general phase approximation in the presence of U(l) gauge- 
pl£I The quantum interference effects can be included by the renormalizatinn of the diffusion constant D and 
othe r par ameters of the theory. In fact it has been argued recently in referencec3 that the phase approximation of 
Eq. ( 3.37 ), along with the effective action of the screened-potential of Eq. ( 3.2S ) can be derived by seeking a temporally 
and spatially varying saddle-point solution Qo{r, t, $) of the action Sq{^) for each <I>(r, t). Quantum interference can 
be treated systematically by considering fluctuations around such saddle point solutions. We note in passing that, 
at criticality, another complication arises: due to the multifractal behavior of the critical eigenstates, the diffusion is 
anomalous, i.e., D becomes dependent on the length and the time scales. These subtleties will be addressed in later 
sections. 

The final step is to substitute Eq. (3.37) into Eq. (3.29) and carry out the functional integral. Taking into account 
the Gaussian fluctuations in $ captured by the effective action in Eq. (3.28), we obtain the central result: 

(G(r,r;r)) « G(r, r; r; 0)6^^^ 

where the Debye- Waller phase-delay factor is 



(3.38) 



Wir) = - 



2 ^ 



(27r)2 



p(q, iLUn)vsciq, ia;„)p(-q, -iw„). 



and Wsc is the dynamically screened interaction implied in the effective action in Eq. (3.28), 

v{q) v{q) 



Wsc(q,iw„ 



1 + v{q)Il{q,uJn) 



In Eq. ( 3.38 ), G(t) corresponds to the SCBA Green's function Ggpir) defined in Eq. ( 3.20 ), 

G(r, r, r) = ^^'^^'o^ ^ e"'""^sign(tj„) = -j/q- 



sin(T7r/ (3) 



(3.39) 



(3.40) 



(3.41) 



After carrying out the sum over Matsubara frequency in Eq. (|3.39| ), the details of which are given in Appendix A, 
we obtain the finite-temperature expression: 



W{t) 



2TTi 



1 



1' 



where 



. ^ 1 f d^q 

= 2 y (2^ 



Dq^ 



Vsciq, -iuj). 



(3.42) 



(3.43) 



The interaction correction to the TDOS is determined by the behavior of the phase- factor W{t), which depends on 
the nature of the dynamically screened interaction Vsc- 
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IV. BULK TUNNELING DENSITY OF STATES IN 2D 



We now derive the TDOS at T = 0. It is necessary to perform the foUowing analytical continuation: 



p{uj) = --lmG{r,r;iujn)\iLOr,^uj+iS = --Im 



dre'^^-^Gir, r, r)e 



(4.1) 



The procedure turns out to be quite non-trivial. Since we could not find discussions of the technique in the literature, 
we elect to include the details of the analytical continuation in Appendix B, where we show, in the limit T ^ 0, 



iy{u!) « -I/O 

TT 



t 



where Uq is the noninteracting TDOS near the Fermi level. Taking the T = limit of Eq. ( 3.42 ), we get 

duj 



W{it) 



2ni 



[f{-^^)-f{^u;)]{l-e-''^% 



(4.2) 



(4.3) 



with the function / given in Eq. (3.43). The term with the oscillatory factor e~*"' averages to zero upon integration 
except for cu ^ 1/t where e~'"* « 1. Therefore we can effectively leave out the e~''^* term and introduce a lower 
cut-off fi/t to the integral: 



W{it) 



1/t 27ri 



(4.4) 



The upper cut-off of the integral in the above equation arises from the fact that the diffusive picture becomes invalid 
at time scales shorter than the elastic scattering time. Using Eq. ( 3.43 ), we obtain, 



W{it) = 



l/ro 



cLu 



It 2^7 (2^) 



Im 



1 



Dq^ — ioj 



(4.5) 



We next turn to the evaluation of the most singular contributions to W and thus to the TDOS for different forms of 
interactions. 



A. Long range Coulomb interaction 



The singularity in the TDOS arises from the ph ysics of dynamical screening. For Coulomb interaction, v(q) 
27re^/g. The dynamical screened interaction in Eq. ( 3.40 ) becomes. 



_kIV 
Dq^-i 



(4.6) 



where n — 2i:e'^dn/ dfi is the inverse screening length at the transition. It is important to note that in the presence 
of disorder, the range of validity for static screening is quite small.ES Since diffusion is a relatively slow process, at 
non-zero frequency, the lo ng d istance singularity associated with the long range Coulomb interaction is not screened, 
as can be seen from Eq. (4.6). In fact, in the region where Dq^ < oj < Dhiq, the effective interaction has the most 
singular form: 



1 



Dnq^ 



Dq < OJ < Dnq, 



which gives the main contribution to the wave-vector integral in Eq. (4.5) 

"i/^" duj 



W{it) 



1 r^^" duj LP 



(4.7) 



(4.8) 



Note that in this region, the diffusion coefficient D is a constant. The anomalous diffusive behavioiS in the regime 
Dq^ 3> does not affect the leading contribution. The remaining frequency integral generates the double-logarithmic 
dependence in time: 
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W{it) 



■ln(-)ln(-), 

To Tl 



(4.9) 



where ri = 1/to{Dk^)'^ , and axx = Ddn/dfj, is the conductivity defined via the Einstein relation. Near the Landau 
level centers, one can show in the SCBA that ti/tq — {\/'^TT^axx){kfaBY /kfi <C 1. This double-logarithmic form is 
the dominant behavior of the Debye- Waller phase factor in the long-time limit. Next order corrections are of the order 
{^/^xxi^l<^tx\^'^{^/'''o)- The contributions fro m all six diffe rent integration regions in the (a;,q)-planc arc discussed 
in detail in Appendix C. Substituting Eq. (19) into Eq. (O), we obtain the zero-temperature TDOS in the Coulomb 
case, 



v{uj) 



2z/o 



t 



1 ln(-!-)ln(-!-) 



(4.10) 



Keeping in mind that in deriving this result we have assumed a frequency independent conductivity axxi i-e. we have 
neglected the quantum interference effects. In general, axx is renormalized by localization effects, of leading 
(l/ca;a;) Ino^To in the unitary ensemble, and by interaction effects of leading order \uujtq in strong magnetic fieldr'^ 
Thus Uxx takes on the frequency-independent SCBA value only if | ln(eTo)| <^ (Txx- 



1. High frequency regime: \ ln(ajTo)| <C <Jx 



In this regime, the weak localization correction to the conductivity can be neglected. If in addition, |ln(ti>To)| <C 
/(Jxx: we can expand the exponential in Eq. ( 4.10| ) to leading order in 1/ axx and obtain. 



■ln(|w|To)ln(|tj|ri) 



(4.11) 



This reproduces the high-field perturbative diagrammatic result of Girvin, Jonson, and Le eO a nd Houghton, Senna, 
and Ying.La For frequencies in the range, y/a^ <C | ln(|ti;|To)| <C axx, the integral in Eq. (4.10) can be evaluated by 
the stationary point/instanton method, leading to a nonperturbative resummation of the double-log divergences in 
Eq. (pll), 



■ ln(|ti;|ro) ln(|cj|ri) 



lyitu) — vq exp 



In the zero magnetic field case, such a nonperturbative resummation of the pert 
carried out by Finkerstein,E3 and recently reexamined using different approaches 
be regarded as an extension of the latter to the strong magnetic field case. 



(4.12) 

double-log divergences was 
Our result of Eq. (4.12) can 



2. Low frequency regime: \ ln(eTo)| 3> <Jx 



However, at 
been shown 
^ Thus, the 



Here the quantum interference effects will, in general, lead to a frequency dependent conductivity, 
the IQHT, the critical conductivity ct^ is finite and of the order of e'^/h. This experimerita,]_facj_ha. 
numerically for both noninteracting electrons and interacting electrons in the HE theory.L2l-L£lL3t!j'E^ 
validity of our analysis, i.e. the structure of the double-log divergence at long-times, can be extended into the regime 
of small u!, provided that axx in Eq. ( 4.10 ) is re placed by the critical conductivity CTc — 0.5/27r. Note that the due to 
the double- log term in the exponent, Eq. ( 1.10| ) implies 



lim W{it) = +00, 



and consequently a zero-bias anomaly in the TDOS 



v{uj = 0) = 0. 



To obtain the limiting behavior of v{lo) for small uj, we expand the sin(wt) factor in Eq. ( 1.1C| ) in a power series in 
Lot. It is important to emphasize that this is possible because of the double-log contribution which makes the time 
integral over e~^^**^ copyerge fast enough such that the TDOS becomes analytic at small uj. The claims made by 
Polyakov and SamokhinES that the TDOS falls off faster than any power law in u; is in fact incorrect. Since the signs 
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of the expansion-coefficients alternate, the series is asymptotic, i.e., it can be infinitely accurate at small uj. To first 
order in w, 



v{uj) = vq\uj\- 



dte 



1 — In(-i-) In(-i-) 

T^a^ ^^Q' ^-^l' 



(4.13) 



Performing this integral, and use the fact that the compressi bili ty is only weakly renormalized, i.e. dn/dii ~ i/q, we 
obtain the 2D quantum Coulomb gap behavior given in Eq. (|l.^) in the introduction, i.e. 



(4.14) 



In contrast to the 2D classical Coulomb gap, the slope Cq is not a universal number, but rather a quantity of quantum 
mechanical origin that depends on the microscopic details of the sample. 



C'c 



1 



1 + *(V2^i^ 



,(27rVe) + (l/87rVc) log^ R 



(4.15) 



where #(x) is the error function and R — \/t\Itq = 1/Dk^tq is a quantity that depends on the degree of disorder. 
The latter can be written in terms of more familiar quantities according to 



R = 



1 (fcyfls)' 



Aw'^ac kfio ' 

where as is the Bohr radius and £o is the zero-field mean free path. 



(4.16) 




FIG. 2. The TDOS in Eq. (1.10) in the case of 



long-range Coulomb interaction, showing the asymptotic 
linear Coulomb pseudogap behavior at small \uj\. The pa- 
rameters are a^x = o^c — O.be^ /h and tq/ti — 10. 



It is easy to verify that the next term in the expansion is of the order uj(ujtq) 



, which is small in this regime. 
The results of numerical int egrat ion of Eq. (|4.10| ) is plotted Fig. 2, which shows t he cro ssover from the high frequency 
behavior described by Eq. (4.12) to the asymptotic linear Coulomb gap of Eq. (4.14) at low frequencies. Since the 
real transition must be governed by an interacting fixed point where the Coulomb interaction strength is finite, we 
conclude that the true asymptotic behavior of the bulk TDOS exhibits the quantum Coulomb gap at the IQHT. 



B. Short-range interactions 



In this subsection, we address the question of how short range interactions, such as screened Coulomb interactions, 
which are irrelevant perturbations at the NIFP in the RG sense, cause depletion of the TDOS near the Fermi 
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level. For simplicity, we focus on the local interactions described by the p rototype short range interacting potential 
v{r — r') = u5{r — r') and v{q) = u. The screened interaction in Eq. ( 3.40 ) becomes, 



1 -I- iiiH—S^— 



(4.17) 



Inserting this expression into the Eq. ( [4.5| ) for W{it), one finds 



Wsriit) 



ZTT 



Im 



1 



1 



Dq"^ — iuj D'q^ — ioj 



(4.18) 



where D' = D + uaxx- In contrast to the long range Coulomb case, the contributions to the g- integral from the 
w > Dq^ regime and the uj < Dq^ regime are now comparable. 



1. Pre-scaling regime 



Let us first ignore the quantum interference effect and focus on the perturbative regime appropriate when | In cotq \ 
y/(Txx- In this case we can treat the diffusion coef ficien t _D as a constant, and the interaction strength u as a marginal 
perturbation (a scale invariant constant) in Eq. (4.18). Carrying out integrations, we arrive at. 



/ dLO-{l 

Jo ^ 



where a is a non-universal dimensionless quantity dependent on the interaction strength, 

1 2 + A 



a = X- 



\7r^axx (1 + A)^ 



(i + IuVi + a) 



(4.19) 



(4.20) 



with A ~ ui/Q. Thus W{it) divergea,in the long time limit only logarithmic ally. T he sit uation is completely analogous 
to the classic X-ray edge problem.cj Taking care of the short-time behavior in Eq. ( 4.19|) as in the X-ray edge problem. 



we obtain, 



Substituting Eq. (4.21) into Eq. 



I^sr(t^) 2 

— — 



W{it) ~= aln(l + i/T-o). 
, we have for the TDOS 

1 



-^^sin(|..) 



t (1 + i/^o)" 



(4.21) 



(4.22) 



is a dimensionless numerical constant. Thus we conclude that a nonuniversal 



where C^, = {2tt) dy—-^j^^ 
power-law suppression of the TDOS prevails in the prescaling regime for short-range interactions. 



2. Scaling regime 

On approaching the scaling regime of the IQHT, it is necessary to take into account the scaling behavior&-pf (1) the 
diffusion coefficient D and (2) the interaction strength u. It is known from the work of Chalker and DanielE^ that the 
multifractality of the critical eigenstates leads to anomalous diffusion in the regime Dq^ > oj. The diffusion constant 
becomes a function of q^ juj: 

D{q,Lo) = D{qL^)-\ (4.23) 



where as before, L,^ = yj D /uj and rj ~ —x^ which is given in Eq. ( pl^ ), and D2 ^ 2 — T] is the multifractal dimension 



The critical conducti vity c Tr. is once again finite and will be taken as scale- independent. The Debye- Waller phase 

I a a' where 

+ log(l + A)] , (4.24) 



factor Wsr(*i) in Eq. ( 4. IE ) is now modified by replacing a a' where 

1 



87rVc 
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with Crj = 1/2 + 2/(4— 3i]). This modification due to the mu ltifra ctal behavior alone would lead to, repeating the 
calculations above, the same behavior of the TDOS as in Eq. ( 4.22| ) in the prescaling regime, except the exponent a 
is replaced by a'. 

Next we r nuat .take into account the fact that u is an irrelevant perturbation. The effective interaction scales to zero 
according tcOllS UcS ^ uuj^+^^ , where —x+ is the scaling dimension for short-range interactions discussed in section 
II and z = 2 is the dynamical exponent at the NIFP. As a result, the quantity a' obeys the following scaling relation. 



The fact that A{u ^ 0,a;) = imphes, together with Eq. ( [4.24D , the leading scahng behavior for a', 
Substituting this result into Eq. ( 4.19| ), we find 



(4.25) 



(4.26) 



(4.27) 



where 7 = x+/z ~ 0.32. That War {it) converges now in the long-time limit should be contrasted to the long range 
Coulomb case and is a consequence of the short range interactions being irrelevant, i.e. 7 > 0. An immediate 
implication is that the TDOS would be finite at zero bias and the level spacing scales as 1/L^ as in the noninteracting 
theory. However, we shall show below that if the bare interaction A is large, it still leads to strong suppression of the 
TDOS at low energies. 



From Eqs. and ( [4.27| ), the TDOS is given by 



dt 



M\^\t) JAXH)(rn/t)-' 



(4.28) 



where i'{0) = i^oe ^^^f < i^q. Performing the integral using the saddle-point/instanton approximation, we find that 
at low frequencies, lo < ujq, the TDOS is given by 



1+ ( ^ 



(4.29) 



where ojq = Tq^{AX/j)~^^'^ is an energy scale. We see that upon scaling, the irrelevance of short-range interactions 
leads to a smearing of the power-law behavior in the perturbative regime, giving rise to a finite and nonuniversal 
TDOS at zero bias. However, although short-range interactions are irrelevant in the RG sense, since iy{Q) ^ vq fo r 
large A, they still lead to strong density of states suppression at low bias. What is remarkable is that Eq. ( 4.29| ) 
predicts an increase of the TDOS with energy that follows a universal power law, with an initial cusp singularity 
for our value of 7. These predictions can, in principle, be tested experimentally by deliberately screening out the 
long-ranged Coulomb interaction using metallic gates or ground planes. 



General interacting potential: v{q) 



It is interesting to consider a general interacting potential of the form v{q) = ^t^- We find that for p < 2, the 
dominant contribution to W{it) in the long time limit comes from the same double-log term as in the case of Coulomb 
interaction (corresponding to p = 1). We conclude that for p < 2, the density of states in the asymptotic w — > 
limit is of the form of the linear gap iy{uj) ^ lo. For p > 2, the phase-delay factor W(i£) approaches a constant 
W^oo > for t ^ tq. In this case the density of states does not vanish, but rather develops a shallow dip at a; = 0, 
where i/(0) = VQe~^°° . The borderline case is that of p = 2, corresponding to the ^-function interaction studied. It 
can be shown from either the p > 2 or the p < 2 side that as p — > 2, a single-log term emerges and dominates the 



contributions in W{it), leading to the power law density of states of Eq. (4.22) in the perturbative regime, and to the 



finite zero bias TDOS we obtained in Eq. (4.29). Details of this analysis can be found in Appendix C. 
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V. TDOS IN QUASI-ID SYSTEMS WITH EDGES 



In the cases studied above, the bulk TDOS does not depend on the Hall conductance. This is in keeping with 
the fact that the bulk diffusion equation is the same with or without time reversal symmetry. The traverse force 
induced by a magnetic field does not affect the diffusive charge spreading. It is well-known that in the rpii-interacting 
theory of the IQHT, the term in the action that depends on a^y is topological and non-perturbative.ll3 However, it 
has been discovered recently that in the presence of edges, the Hall conductance enters measurable quantities even in 
the perturbative limit. The topological term gives rise to a tilted boundary condition for diffusion and at mittce subtle 
levels affects the quantum interference processes. For example it was shown by Khmel'nitskii and Yosefin,Eil and by 
Xiong, Read and StoneE3 that the mesoscopic conductance fluctuations in phase-coherent samples become dependent 
on the Hall conductance in the presence of edges. More recently Shytov, Levitov, and Halperin demonstrated that the 
I-V curves for edge tunneling into the ID Luttinger liquid like edge excitations of fractional quantum Hall liquids can 
be obtained from the point ojLview of bulk composite fermions by using a similar phase approximation in treating the 
effects of gauge fluctuations. e3 There, the a^y dependence in the exponent of the power-law tunneling conductance 
also arises from a boundary condition of the source current at the edge of a semi-infinite sample. Our case differs 
from and is simpler than that of the composite fermions in the sense that we need to consider only pure potential 
fluctuations mediated interactions in the integer quantum Hall regime. 

To study how-the physical boundaries bring the Hall ratio into the bulk TDOS, we consider, instead of the half 
plane geometry,Ea a quasi-lD sample with its length L much greater than its width W, exposing two reflecting edges 
along its width. This condition can be-*ealized experimentally in the long Hall bar geometry, and is the same as that 
considered by Xiong, Read and StoneE3 in their study of the edge effects on mesoscopic conductance fluctuations in 
strong magnetic fields. 

Because the incident current is at an angle with the reflecting edges, the presence of the magnetic field affects the 
diffusion process through a modified boundary condition which depends on the Hall ratio — cFxy/o'xxj 



[dn + iHdt] p = 0. 



(5.1) 



Here the subscript n denotes the directions normal to the edge and t denotes the tangential direction. Strictly 
speaking, if the new boundary conditions (5.1) are taken into account, the diffusion propagator as well as the screened 
interactions will depend on the Hall conductivity. Such effect is minimal if the sample is wider than it is long and 
becomes pronounced only in the quasi-lD limit when L ^ W . For simplicity, we consider the case of the 5-function 
interaction encountered in the previous section and the geometry of an infinite strip with hard walls at j/ = and 
y — W. We also limit ourselves to the perturbative regime, and neglect scaling corrections to the conductivities and 

the interaction strength. 

The diffusive modes |tiiat are solutions of the diffusion equation (3.34) and satisfy the boundary condition ( ^.l| ) can 
be obtained as follows£3 



-cos( — )-»7«gsin(— ) 



f>L ^ ^^^^qx-^^„qy ^ for fc = 0. 



for fc 7^ 



(5.2) 



Since the boundary condition is not self-adjoint, there are also a set of right eigen-functions (with the same eigen- 
values) that satisfy the boundary condition under parity transformation: 



[dn - iHdt] dp"' = 0. 



(5.3) 



They are given by, 



— n, p^l^ 



kn ,kTTy. . .kiry 

— cos( ) -|- i^HQ sm( ) 

W ^ W W 



ao 



for fc 7^ 



(5.4) 



The {(f) } are the eigenmodes of the Laplacian operator. 



(5.5) 



where 
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A, 



q.k 



D 



2^2 



, for fc 7^ 
for fc = 0. 



(5.6) 



Using the bi-orthogonality relation and the completeness condition, we can express the Debye- Waller phase factor 
W{it) in Eq. (^|^) in terms of the sum over the eigenmodes in the transverse channels: 



W{it) = i 



— -^0^,,(r)[</.^,(r)]*2Re. 

1/* y^""' q,k 



1 



(5.7) 



In the limit Dt ^ W"^ , contributions from A; > modes can be ignored. For the consideration of i^(w), this condition 
translates into W, i.e. hco » D/W'^. Strictly speaking, W{it) becomes dependent on the spatial position, but 

we do not expect any spatial singularity. It is therefore justifiable to average W{it) over the entire strip. We obtain. 



Wiit) 



37r2 hy/D^ 1 + A 



(5.8) 



where Di ~ (1 + -ffj)D. Thus , th e phase factor is dominated in quasi-lD by the ^/t divergence in the long-time limit. 
Rewrite the prefactor in Eq. (5.S) as 



/ujb = 



(5.9) 



we have W{it) 
by, 



/LUst — y^ujETQ. Upon substitution of W{it) into Eq. (O), the TDOS in the quasi-lD case is given 



dt 



Approximating the time-integral by the stationary point /instanton solution, we find that for a; <C w^. 



(5.10) 



(5.11) 



Thus, the TDOS is strongly suppressed at low energies in quasi-lD systems in a manner that is se nsitive to the applied 
magnetic field through the energy scale ujb- If one naively extends the behavior of Eq. (5.11) to arbitrarily small 
frequencies, one would have concluded that the TDOS goes to zero on the Fermi surface faster than any power-law, 
if the renormalization of th e conductivity and the interaction strength at low energies is ignored. This is in fact 



incorrect, because Eq. (5.11) is only valid at intermediate frequencies. Since Debye-Waller factor e 



-W{tt) 



converges 



fast in the long time limit, the low energy behavior of v^id{uj) is actually described by an asymptotic series expansion 
in powers of lo. Wc find. 



i^B ^} (2n 

n— 



lub 



The TDOS is therefore dominated by the linear term near zero bias, 

t'qlD(w) ~ fO— = S\uj\- 
UJB 

The magnetic field dependent slope is given by 



9n^ 



1 



16m2 p^, 



-(l+A)^ 



(5.12) 



(5.13) 



(5.14) 



where p^x — o'xxi^ + 7h) is the dissipative resistivity. Note that this result is valid at small frequencies such that 
\uj\/ujb > e'^^/'^l'^l, i.e. for \uj\ < 0.12^^. In Fig. 3, we plot the TDOS obtained by numerical integration of Eq. ( |5.10| ) 
as a function of lo. The asymptotic linear pseudogap behavior of the TDOS at low bias is shown in the inset. 
Therefore, we conclude that the TDOS of a quasi-lD quantum Hall strip with reflecting edges exhibits a linearly 
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vanishing pseudogap near the Fermi level, with a slope proportional to in the perturbative regime. Within the 
SGBjA, the values of both Gxx and Uxy at the center of the Landau levels are proportional to the Landau level index 
A^.Ea This leads to a Hall ratio of order one and a longitudinal resistivity pxx ^ ^/N ^ B. We see that for a fixed 
interaction strength u, the slope is inversely proportional to the magnetic field s ^ 1/B. ft is interesting to remack 
that bulk tunneling measurements under the quantum Hall conditions using time-domain capacitance spectroscopyEil 
indeed reveals a linearly vanishing pseudogap TDOS with a slope that scales with 1/ B. In addition, the measure slope 
of the tunneling pseudogap oscillates weakly as a function of filling fraction, which mimics the oscillatory behavior of 
the SCBA conductance. However, although the experimental setuptj allows screening of the Coulomb interaction by 
the metallic gates (electrodes) such that the interactions may be short-range, it is not clear at the present if the the 
sample used can be qualified effectively as being quasi-one-dimensional. It is also interesting to note that the results for 
the quasi-lD bulk TDOS in Eqs. (5.13) and (5.11) depends strongly on the interaction strength u. This is in contrast 
to the case of tunneling into a single fractional quantum Hall edge in the composite fermion description, wiic;re the 
TDOS was found to be a power-law with an exponent that depends only weakly on the interaction strength 




FIG. 3. The TDOS in quasi-lD strips with edges ob- 
tained from Eq. (5.10) for short-range interactions. The 



inset shows the asymptoti c low energy behavior of a linear 
pseudogap given by Eq. (5.13). 



VI. DISCUSSIONS: 



QUANTUM COULOMB GAP AND DYNAMICAL SCALING OF TRANSITION 

WIDTH 



Understanding how interactions introduce new physics in the single-particle sector near the IQHT is an important 
step towards a more complete understanding of the interplay between disorder and correlation and its effects on the 
transport properties in the quantum Hall regime. A great part of this article is devoted to demonstrating how various 
types of interaction-induced anomalies in the TDOS at low energy are likely to occur so long as the 2D conductivity 
is finite which happens near the quantum Hall transitions. Our basic finding is that, in the presence of disorder, 
the range of validity for static screening of the Coulomb interaction is very small, whereas at finite frequencies the 
diffusive dynamics is to slow to effectively screen out the Coulomb interaction at long distances. This leads to a 
Debye- Waller phase-delay factor in the averaged single-particle Green's function that vanishes in the long-time limit. 
As a result, the TDOS in the extended regime comes to resemble that in the localized regime, exhibiting a linearly 
vanishing quantum Coulomb gap behavior. 

It is important to emphasize the coexistence of the va nishin g Coulomb gap in the tunneling DOS with a finite 
thermodynamic DOS. The double-log singularity in Eq. (4.10), arising from the correlations of the single particle 
phases, will not show up in the two-particle density-density correlation function that determines the compressibility 
or the thermodynamic DOS iri-ihe static limit. This point has recently been demonstrated explicitly in the self- 
consistent Hartree-Fock theory,EZI where it is shown that the charge redistribution induced by a test charge inserted 
into the quantum Hall critical state is consistent with the presence of a finite screening length. The finite critical 
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conductivity then imphes that the uniform diffusion constant must be finite. 

It is also important to understand how the depletion of TDOS relates to the larger issue of dynamical scaling near 
the IQHT. This is more challenging due to the possible existence of several different time scales that control different 
dynamical processes: charging, charge spreading, and inelastic phase breaking, etc. While a linearly vanishing single- 
particle DOS in 2D does imply, through the quasiparticle level spacing A ^ a dynamical scaling exponent z — 1, 
it has not been shown how this z, which obviously controls the charging dynamics, is also the one that enters the 
conductivity scaling at finite temperatures or frequencies in the transport measurements. 

The conventional wisdom, at least for d.c. transport, has been that the hopping transport in insulators is determined 
by the single-particle DOS, whereas the diffusive transport in metals is directly related to the thermodynamic DOS 
(compressibility). Generalized to finite frequencies, this would imply that the dynamical aspects of the transport in 
insulators are controlled by that of charging, and in metals by that of dif&isive charge spreading. In ordinary disordered 
metals, this problem can be better quantified. As noted by Finkel'steinES in the calculation of the conductivity, there 
is a great degree of cancelation between the corrections to the Green's functions (which causes the anomalous behavior 
in the TDOS) and to the vertices. Therefore the interaction effects that cause the depletion of the the single-particle 
DOS h'ii-u) do not appear to influence directly the d.c. transport properties. At finite frequencies, the coupled scaling 
equations in the RG calculation based on the nonlinear-cr model involve the conductivity, the interaction strength, 
and the frequency renormalization that enters in the diffusion kernel as l/{Dq- — iZ^oj). Note that Z^^ = 1 
corresponds to the non- interacting case and a dynamical scaling exponent z^j — dim[-D] -I- 2 = d in d-dimensions. In 
conventional Anderson-Mott metal-insulator transitions in d = 2 -I- e dimensions ,lZI the critical conductivity is zero. 
The lack of quantum diffusion at the transition is accompanied by the interaction induced frequency renormalization, 
i.e. Z^ ^ . For both the spin-scattering as well as the spin-polarized case, one-loop calculations give C, = e/2. 
As a result, the dynamical scaling exponent relevant for metallic transport departs from the noninteracting value, 
z^ = d — ^. It is very difficult to extend the same quantitative analysis to the quantum Hall problem because the 
perturbative approach is no longer valid due to the presence of the topological 0-term in the nonlinear-tr model action. 
However, the fact that both dn/dfi and the conductivity are finite at the quantum Hall transition in 2D ensures that 
the frequency in the diffusion propagator remains unrenormalized and its associated exponent z^j = 2. 

The above analysis conveys a simple but important point, i.e. since it is and not the single-particle DOS 
that enter s sca ling and controls the dynamics of the diffusive transport from the metallic side, within the existing 
frameworkc3E3, it is natural to suspect that the suppression of the TDOS plays no role in the dynamical scaling 
behavior of the conductivity. Therefore, the linear Coulomb gap may appear not by itself an explanation as to 
why z = 1 at the quantum Hall transitions. In the following, we argue that it is indeed the interplay between the 
quasiparticle inelastic dephasing (level broadening) and the level spacing that controls the transition width, contrary 
to common perceptions. 

Since the relevant phenomenon here is transport, this state of affairs naturally translates into the question of which 
exponent determines the dephasing length L^p. In the standard procedure,E3 the system is divided into L^p x phase 
coherent blocks. Transport within each block can be described by phase coherent transport from the underlying 
noninteracting theory and the relevant conductivity is given by the disorder-average over the phase-coherent blocks. 
The outcome is that the scaling variable for the conductivity becomes L^/i^ in the presence of interactions instead of 
L/^ in the noninteracting case, where L is the sample size and ^ is the localization length, 

= x'^ (t) ■ 

The scaling function F{x) has the limiting behavior. 

The conducting critical regime at L^p <^ ^ and the insulating regime at ^ ^ are separated by a crossover at ~ ^, 

where the scaling variable in Eq. (6.2) is of order one, giving rise to a transition width S* ^ Lp^^^^"". Physically, 
the transition width can be viewed as the width of the energy window of states whose localization length exceeds the 
phase coherence length. In the language of quantum critical phenomenon, the latter corresponds to the width of the 
quantum critical region. This is a generic property associated with the quantum critical point. The only peculiarity 
is that on either side o£,the quantum Hall critical point the ground state are insulators which, drawing analogy to 
quantum spin systems J23 are quantum disordered (see Fig. 4). The renormalized classical regimes (metallic phases) 
are absent. nnn 

Now, we examine the conventional view of finding L^. For a generic quantum phase transitionQQa the critical 
regime is characterized by the only time scale h/T, and thus the dephasing time ^ XjT . If quantum diffusion is all 
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that matters, the associated length scale, i.e. the thermal diffusion length would be set by Lt = (i^r^)^/^ 
leading to a thermal exponent zt = 2. Similarly, the length associated with a finite frequency is — {Dfiijj)~^^^ such 
that the dynamic exponent z,^ — 2. The conventional approach is to identify Lt,u with the quasiparticle dephasing 
length in the scaling function in Eq. ( |6.lD . As a result, the transition width in this picture is determined by the 
crossover of length scales set by Lt,ui ^ indicated in Fig. 4 by the dashed line, which leads to 

Cuid-be-(r^/^"^'-,^'/^'^'^'-)- (6.3) 

As emphasized in the introduction, with the values v ~ 2.3 and zt = z^^ ^J2~ Eq. (|6.3[ ) does not agree with the 
scaling behavior of the transition width measured by transport experiments J13'll3 From the theoretical point of view, 
'^wouid-bo would be the width of the phase coherent, diffusive metallic transport regime, provided that Coulomb 
interaction effects are not too strong to induce single-particle localization of the quasiparticle states. 



T,CO 



Quantum Critical Region 




FIG. 4. Schematic phase diagram near the quan- 
tum critical point {Ec) of the IQHT. The dashed line 
describes the "would-be" crossover between phase coher- 
ence, diffusive metallic transport of the quantum critical 
region and the localized quantum disordered regime at 
Lt,i^ ~ C- The associated thermal and frequency expo- 
nents are zt = z^j = 2. The solid line describes the true 
crossover from diffusive, metallic to single-particle, insu- 
lator-like transport (shaded area) that talces place when 
Lip ~ The observed transition width, 5*, is narrower 
than ^wouid-bc has a scaling behavior controlled by 
the dynamical exponent z = 1, arising from the quantum 
Coulomb gap behavior. 



We now point out the problem with this picture which is commonly used to describe metal-insulator transitions. 
The use of Lt.u) as in Eq. (6.1) completely ignores the important Mott physics in the single-particle sector, namely 
the tendency toward the single-particle insulator behavior induced by Coulomb interaction. Physically, the inelastic 
dephasing time is related to the interaction induced quasiparticle decay rate or level broadening F h/Tip. The 
dephasing length L^p^ on the other hand, can be determined only if the underlying transport mechanism is specified, 
ballistic or diffusive, insulating or metallic. Clearly, diffusive metallic transport is possible only if there is a significant 
overlapping of the quasiparticle levels within F, i.e. the level broadening must be larger than the interaction induced 
level spacing, F > A. In the opposite limit, F < A, the single-particle states are essentially gaped as a result of 
electron-electron interaction and the transport would be controlled by the localization in the single-particle sector 
similar to that in a Mott-insulator. Diffusion would be impossible in this case and one expects variable range hopping 
to be the dominant mechanism of transport. It is therefore more appropriate to define the quasiparticle dephasing 
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length by the condition T ~ A{Lip). The presence of the quantum Coulomb gap, Eq. (4.14), in the critical regime 
implies that near the Fermi level, 

A(i^) ^ -j^, (6.4) 



which leads to ^ ^ 1/T and the dynamical exponent z = 1. Note that L^/Lt ^ T^^/^, i.e. L^p ^ Lt at low 
temperatures. However, it is Lip that controls the crossover to the quantum disordered insulating regime since when 
Lip > ^, the level spacing within a ^ x ^ volume becomes larger than the level broadening and the diffusive metallic 
transport becomes impossible. It must be replaced by hopping dominated transport similar to the|-Goulomb blockade 
regime in quaiitum dots. This part of the physics has been emphasized by Polyakov and ShklovskiiE^ and by Polyakov 
and SamokhinEa in terms of the classical Co ulom b gap. 

Using this L^, the scaling function in Eq. ( |6.l| ) now describes the true crossover line, set by ~ ^ (the solid line in 
Fig. 4), that separates the diffusive, metallic transport from the single-particle, insulator-like transport. The scaling 
behavior of the width of the critical conducting regime is therefore given by 

5* {T,Ljy (6.5) 

with z = 1 as observed in transport measurements. In this theory, the presence of the quantum Coulomb gap behavior 
is central to the emergence of the z = 1 scaling of the transition width. It eliminates the difhculty associatejd with 
invoking the classical Coulomb gap or the bare charging energy due to the unscreened Coulomb interactiorH'ES which 
is only valid deep in the insulating regime. 

As shown in Fig. 4, the single-particle DOS obeys the quantum Coulomb gap behavior in the quantum critical 
regime, whereas deep in the quantum disordered, insulatiua regime, it is expected that the Coulomb interaction 
reinstates the classical Coulomb gap of Efros and ShklovskiiE3 Thus it is remarkable that in the presence of Coulomb 
interactions, the crossover between the quantum critical and the quantum disordered regimes is accompanied by a 
cros sover in the behavior of the TDOS — from the quantum to the classical Coulomb gap. Comparing Eqs. (|l.4| ) and 



(|1.8D, we see that the crossover is simply described by a crossover in the slope of the linear gap near the Fermi level. 
Such a crossover should, in principle, be detectable experimentally by sitting at a fixed distance to the critical point 
of the transition. In this case, as the temperature or frequency is lowered, one should observe a linearly vanishing gap 
with an initial nonuniversal slope that turns into a universal number in the low temperature/frequency limit. 

We emphasize that the linear quantum Coulomb gap behavior results from the combined effects of (i) two- 
dimensionality, (ii) long-range Coulomb potential, and (iii) quantum diffusion, i.e. a finite conductivity at T = 0. It is 
expected to pertain to other metal-insulator transitions in 2D amorphous electron systems, provided that the critical 
conductivity is finite. The physics discussed here is quite generic of thcj-2|D disordered metal-insulator quantum critical 
point. A recent example is the 2D B = metal-insulator transition.Ej In that case. Fig. 4 needs to be modified to 
include the renormalized classical, i.e. the metallic region. It is our hope that the present work will stimulate further 
experimental investigations on the nature of dynamical scaling in the quantum Hall effect and in other metal-insulator 
transitions. 
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APPENDIX A: THE MATSUBARA SUM 



In this appendix, we carry out the discrete frequency sum in the phase delay given in Eq. (3.39): 

Wir) = ^J2^1-e^-"^]f{M), (Al) 

n 

where w„ = 2TTn/f3, n — 0, ±1, ±2 ... is the boson Matsubara frequency and 
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/(I^nl) = 



1 + v{q)iyQ-. 



To perform the Matsubara sum, we first separate the positive and the negative frequencies by writing 

W{t) = W+{t) + W-{t), 

where 



ri>0 



«<0 ' 



]/(-^n). 



Next, we define a function on the complex plane, 



1 - e^^ 
e^^ - 1 



(A2) 

(A3) 

(A4) 
(A5) 

(A6) 



We consider the integrals of F{z) along contours in the upper (Ci + C2) and the lower (C3 + C4) half-plane as shown 
in Fig. 5. The results are given by, respectively: 



F+{z)dz ^ ^ e"^"^]/(u;„) + 27ri x residues from f{-iz) 



C1+C2 



/3 

2m 



Ji>0 



C3+C4 



F+(z)dz = — y [1 - e*""^]/(-tJ„) + 27ri x residues from fUz) 
p ^ — ' 



(A7) 



n<0 



It can be shown straightforwardly that the integrals of F{z) along both semi-circles C2 and C4 (with \z\ — R) vanish 
at infinite radius R — > 00, provided that < r < /?. Since residues of f{—iz) lie in the lower half-plane while those of 
f{iz) lie in the upper half plane, they do not contribute to the contour integrals as we defined. Therefore, summing 
up the integrals along Ci and C3, we obtain, 



W{t) 



+00 



de 



1-e^ 



s/'^ - 1 2t 



:[fi-^e)-f{^e)]. 



(A8) 



Note that since the fermionic Green's functions are anti-periodic in r, i.e. G{t + f3) = — G{t ), it implies that, through 
Eq. ( |3.38 ), a periodic phase factor W{(3 + t) = W{t) which is indeed satisfied by Eq. ([A^). 
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APPENDIX B: ANALYTICAL CONTINUATION 



In this appendix, we describe one of the technical subtleties encountered when taking the analytical continuation in 
Eq. (O). We show how to obtain G{ui + iO~^) directly from the time ordered imaginary-time Green's function G{t) 
by analytical continuing r — > it + 0"*". We begin with the Fourier transform of the fermion Green's function, 



/■/3 

G{iuJn) 



f dTe"^"^G{T), (Bl) 
Jq 



which satisfies the anti-periodic boundary condition G{t + /3) = —G{t). As a result G{ioJn) is non-zero only for odd 
Matsubara frequencies, i.e. for ujn = i^Iiiii^. 

To perform the integration in Eq. (Bl) and the analytical continuation to real frequency, we extend r to the complex 
z-plane with Re[z] = t and ilm[2:] = it. We seek to analytically continue the integral in the segment bounded by (0, /?) 
on the real axis to integrals along the vertical axis at r = 0, (3. To this end, consider the closed-path integral along 
the contour shown in Fig. 6, chosen to lie in the upper half-place for cun > 0. Since G(r) is non-analytic at r = 0,/?, 
the vertical segments of the contours are shifted infinitesimally such that < Rez < /3. The analytical continuation 
is possible when G{z) is analytic and has no poles encircled by the contour. 



dze*"""G(z) = 0, (B2) 
Since the integral along the j^j ^ oo segment of the contour in Fig. 6 vanishes for Un > 0, we have 

/•oo 

G{ioJn)=i die-'^"*[G(it + 0+)-G(it + /?)]. (B3) 

^0 

Using the anti-periodic property, we obtain 

roo 

G{iuj„) = 2i dte-l""l*G(it-|-0+), (B4) 
Jo 

where we have included the result for Un < 0, in which case, the integration contour was chosen to lie in the low 
half-plane. 

Next we take the analytical continuation in frequency: itUn —^Lu + iO"*" and obtain, 

/>00 /"OO 

G{uj + iS) dte"^'G{it). ~2i e*"*Goe-^('*\ (B5) 



where we have used Eq. (3.38) for the Green's function G(it) i n our semi-classical phase approximation. Substituting 
the expression of the SCBA Green's function Go in Eq. ( ^3.41 ), we obtain the TDOS at finite temperatures, 

= -llmG(. + ^S) = .o| f J^^^.-^i-l (B6) 
TT (j J^^ smh(7rVp) 

Note that an overall factor stemming from the Fermi distribution function has not been included in the definition of 
I'iuj) at finite temperatures, since it is, at any rate, unimportant at low temperatures. 
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FIG. 6. The contour of integration used in Eq. ( p32| ) 
for carrying out the analytical continuation. The real axis 
corresponds to the imaginary time r, whereas the imagi- 
nary axis is labeled by the real time t. 



APPENDIX C: SIX INTEGRATION REGIONS 



In this appendix, we perform the integrations over the frequency and wave- vector in Eq. (4.5) to obtain the phase- 
elay factor W{it). We consider the interaction pc 
v{r) = u/rP in real space. We can rewrite Eq. (4.5 



delay factor W{it). We consider the interaction potential of the general form v{q) = u^^rj in momentum space or 



as 



W{it) 



dw f dq 



i/t 



27ri / 47r2 



Q('?,w), 



where the integrand with p = i/qu is given by 



2uuj{2Dq^ + DqPK 



[{Dq^f 



(Cl) 



(C2) 



There are six characteristic integration regions differed by the ordering of Dq^, Dq^n^^P and uj. The details of the 
integration are discussed below. We find for all p < 2 the leading contribution is the double-log term appeared in the 
long-range Coulomb case (p=l) in Section IV. A. For p > 2 the leading contribution is of the form Woo ^ with 
(5 > 0, and (5 ^ as p — > 2. 



1. The regime Dq'^ < to < Dq^n'^-P 

The condition Dq^ < uj < Dq^n^^P requires uj < Dk^ ior p < 2; and uj > Dk^ for p > 2. In this case, the integrand 
Q in Eq. ( |C2|) can be approximated by, 

Q(p,c.)«— (C3) 

We get for p < 2, 

—Q{q,Lu)^~ — - logwT, C4 

47r^ 27r Dvo p 

where we have introduced r, via 1/ts — Dk^. Performing the remaining uj integral, we obtain 
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where = tq/Dk"^. 

For > 2, due to the requirement that oj > Dk^, this regime does not have any time-dependent contribution in the 
long-time limit when t > Ts- 



2. The regime Dq^ < Dq^K^'^ < uj 

This regime requires q < qo{io) = min[K, (wTs)^/^^] for p < 2 and k < q < {uiTsY/^k, for p > 2. The latter 

case requires w > l/rg. This means that for p > 2, this regime does not contribute in the long-time limit in a 
time-dependent way. 



For j3 < 2, we have 



such that 



/ 



Q{q,uj)^\vou''Dq^P-\ (C6) 



^Q{qM = ^y.u'D / dqq'^-' = \u'd'^^. (C7) 



Usually t > To > Ts, therefore co < ^. We get 

C ^ J S^^''"^ = i-oU^i^^K^^r^og(i/ro) = i^r.«nog(^/ro). (C8) 

This single-log term is subleading when compared to the leading double-log contribution in regime 1 in the long-time 
limit. 



In this case, we have 



3. The regime Dq^k'^'P < Dq^ < uj 



(C9) 



The limits for the q,LJ integrals are k < q < \/ij/D, u > Dk^ = l/r^ for p < 2; and 5 < min[/t. ^^uj/D] for p > 2. 

The p < 2 case is of no interest in this regime since the lower cut-off of the frequency integral is l/r^ and is 
time-independent in the limit t > Ts- 

We now (liscniss the p > 2 case. For most o f the physical systems the mean-free path £ is greater than the screening 
length K, therefore min{K, y/iv/D} = y/iv/D. The integrals can be carried out according to, 



d^q 2Du 1 
— Q{q,oj) = 



1/^" duj 



l/f 



TT W 

U 



2u 

qdqtf = 



1 



tt^DpI"^ p + 2p-2 



TT p + 2 DP/"^ 

1 1 



(to) 



p/2-2 



_P/2-l 
.^0 



(CIO) 
(Cll) 



4. The regime Dq^K^'P < uj < Dq^ 

In this regime, we have 

^i^u 1 
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The limits of the integrations are: (w/D)^/^ < Q < KiiOTsY^'''; uj > 1/ts for p < 2 and uj < 1/ts for p > 2. The 
p < 2 case only produces time- independent contributions in the limit t > t^. 
For the case of p > 2, it is straightforward to obtain 



2uuj I ^ ^' 1 2u 1 

1/2 P 



-—^Q{q,uj, - „3 , 



1 



1 



duj r dPq , , u 1 
2n J An^ tt^ — p 



pTgK' 



P 



(3p-6)(«;^)6(t,/to)6/p 



[l-{To/tf-'/^] 



(C12) 
(C13) 



where ^ = (2Dto)^/2 is the mean free path. 



5. The regime uj < Dq^K, ^ < Dq 



In this regime, the integrand Q in Eq. ( |C2| ) can be approximated by 



(C14) 



For w < I/tq < 1/ts, the limits of the integrations are q > k for p < 2 and [loTsY/pk < q < k for p > 2. 
In the case of p < 2, we have 



fd^q 2UCU p 1 



2UL0 1 1 
D^TT 6-p K^-P 

ut^kP 



[''^°dw [ d'q , _ Ul 

2^ J 4^^^'^'")-27rS2(6-p) 

On the other hand, for p > 2, we have 
leading to 



[l-(roAr 



2uiLj 1 
D^n (6 



(WT5)(6-P)/P 



i/""" duj r d^q u 1 pt,.kP 

2^7 i^^^'^''^^ ~ ^ 6 - p (3p - 6)(K^)6(r«/To)6/p 



i/t 



[l-(ro/t)3-6/P] 



r2 K.P 



2^Vo^(6-p) 



[l-(roAr 



(C15) 
(C16) 

(C17) 



(C18) 



Finally, in regime six, we have 



6. The regime uj < Dq < Dq^K 



2uj 



vn{Dq 



2\Z ■ 



The limits of the integrations are {uj/DY/'^q < k for p < 2 and q > n iox p > 2. 
In the p <2 case, integrals give 



<fq 
in 

duj f d^q 

i/t 



1 UJ 

2nJ 4.2^^'^'")^ 8^ 



logit/ro) 



2D^. 



-Al/rl-l/t^ 



(C19) 



(C20) 
(C21) 
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The p > 2 case gives, on the other hand, 



1/t 



27r / 47r2 



1 



167r2j/oD 



(C22) 



Note that the (5-potential considered in Section IV. B corresponds to the p ~ case. The result in Eq. (4.21) can 
be obtained from either the p > 2 case or the p < 2 case by taking the limit p ^ 2 using lim^^g ^{1 — y^) = —Iny 
and limp^2(j5 — 2) In(ri) = In(i'ou). 
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